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Abstract 



■ We present a detailed analysis of multiple production of the lightest CP-even 

> ; 

Higgs boson (h) of the Minimal Supersymmetric Standard Model (MSSM) at hign- 
in , _ 

fT) ■ energy e + e colliders. We consider the production of the heavier CP-even Higgs 

: 

| boson (H) via Higgs-strahlung e + e — > ZH, in association with the CP-odd Higgs 

00 

Q^ ■ boson (A) in e + e — > AH, or via the fusion mechanism e + e — » v e v e H, with H 

subsequently decaying through H — > hh, thereby resulting in a pair of lighter Higgs 
O 1 bosons (h) in the final state. These processes can enable one to measure the trilinear 

d : 

Higgs couplings Xnhh and Xhhh, which can be used to theoretically reconstruct the 



Higgs potential. We delineate the regions of the MSSM parameter space in which 
these trilinear Higgs couplings could be measured at a future e + e~ collider. In our 
calculations, we include in detail the radiative corrections to the Higgs sector of the 
MSSM, especially the mixing in the squark sector. 
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1 Introduction 



The Higgs potential of the Standard Model (SM), which is crucial in implementing the 
mechanism of spontaneous symmetry breaking, contains the unknown quartic coupling of 
the Higgs field. As a consequence, the mass of the only Higgs boson in the SM, which is 
determined by this quartic coupling, is not known pQ . If a Higgs boson is discovered and 
its mass measured, the Higgs potential of the Standard Model can be uniquely determined. 

On the other hand, supersymmetry is at present the only known framework in which 
the Higgs sector of the Standard Model (SM), so crucial for its internal consistency, is nat- 
ural |2J. The minimal version of the Supersymmetric Standard Model (MSSM) contains 
two Higgs doublets (Hi,H 2 ) with opposite hypercharges: Y(H\) = —1, Y(H 2 ) = +1, so 
as to generate masses for up- and down- type quarks (and leptons), and to cancel gauge 
anomalies. After spontaneous symmetry breaking induced by the neutral components of 
Hi and H 2 obtaining vacuum expectation values, (Hi) = vi, (H 2 ) = v 2 , tan/5 = v 2 /vi, 



the MSSM contains two neutral CP-eve 
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(h, H), one neutral CP-odd (A), and two 



charged (iP^) Higgs bosons [1 . Although gauge invariance and supersymmetry fix the 
quartic couplings of the Higgs bosons in the MSSM in terms of SU{2)l and P(l)y gauge 
couplings, g and g', respectively, there still remain two independent parameters that de- 
scribe the Higgs sector of the MSSM. These are usually chosen to be tan/? and m^, the 
mass of the CP-odd Higgs boson. All the Higgs masses and the Higgs couplings in the 
MSSM can be described (at tree level) in terms of these two parameters. 

In particular, all the trilinear self-couplings of the physical Higgs particles can be pre- 
dicted theoretically (at the tree level) in terms of tua and tan/3. Once a light Higgs boson 
is discovered, the measurement of these trilinear couplings can be used to reconstruct 
the Higgs potential of the MSSM. This will go a long way toward establishing the Higgs 
mechanism as the basic mechanism of spontaneous symmetry breaking in gauge theories. 
Although the measurement of all the Higgs couplings in the MSSM is a difficult task, 
preliminary theoretical investigations by Plehn, Spira and Zerwas [3], and by Djouadi, 

1 Whcn unambiguous, we shall denote the CP-even Higgs particles as h and H. 
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Haber and Zerwas [4J (referred to as 'DHZ' in the following), of the measurement of 
these couplings at the LHC and at a high-energy e + e~ linear collider, respectively, are 
encouraging. 

In this paper we consider in detail the question of possible measurements of the trilinear 
Higgs couplings of the MSSM at a high-energy e + e~ linear collider. We assume that such 
a facility will operate at an energy of 500 GeV with an integrated luminosity per year of 
£ int = 500 fb" 1 !5J. (This is a factor of 10 more than the earlier estimate.) In a later 
phase one may envisage an upgrade to an energy of 1.5 TeV. Since the 'interesting' cross 
sections fall off like 1/P 2 , the luminosity should increase by a corresponding factor. An 
earlier estimated luminosity of 500 ftT 1 at 1.5 TeV may turn out to be too conservative. 

The trilinear Higgs couplings that are of interest are Xphhi ^hhh, and Xhaa, involving 
both the CP-even (h, H) and CP-odd (A) Higgs bosonso The couplings \uhh and X hhh 
are rather small with respect to the corresponding trilinear coupling in the SM (for 
a given mass of the lightest Higgs boson nth), unless is close to the upper value 
(decoupling limit). The coupling XhAA remains small for all parameters. 

Throughout, we include one- loop radiative corrections jB] to the Higgs sector in the 
effective potential approximation. In particular, we take into account the parameters A 
and /j,, the soft supersymmetry breaking trilinear parameter and the bilinear Higgs(ino) 
parameter in the superpotential, respectively, and as a consequence the left-right mixing 
in the squark sector, in our calculations. We thus include all the relevant parameters of 
the MSSM in our study, which is more detailed than the preliminary one of DHZ. 

For a given value of m^, the values of these couplings significantly depend on the soft 
supersymmetry-breaking trilinear parameter A, as well as on fi, and thus on the resulting 
mixing in the squark sector. Since the trilinear couplings tend to be small, and depend 
on several parameters, their effects are somewhat difficult to estimate. 

The plan of the paper is as follows. In Section 2 we review the Higgs sector of the 

2 These are not the only couplings that occur in the Higgs potential. However, these are the only ones 
which could possibly be measured at future colliders. 



3 



MSSM, including the radiative corrections to the masses. The trilinear couplings are 
presented in Section 3. In Section 4 we review the possible production mechanisms for 
the multiple production of Higgs bosons through which the trilinear Higgs couplings can 
be measured at an e + e~ linear collider. In Section 5 we consider the dominant source of 
the multiple production of the Higgs (h) boson through Higgs-strahlung of H, and through 
production of H in association with the CP-odd Higgs boson (A), and the background to 
these processes. This source of multiple production can be used to extract the trilinear 
Higgs coupling X Hhh . 

Section 6 deals with a detailed calculation of the cross section for the double Higgs- 
strahlung process e + e~ — > Zhh. This process involves the trilinear couplings Xnhh 
and Xhhh of the CP-even Higgs bosons (h, H). In Section 7 we consider the differ- 
ent fusion mechanisms for multiple h production, especially the non-resonant process 
e + e~ — > fevjih, for which we present a detailed calculation of the cross section in the 'ef- 
fective WW approximation'. This process also involves the two trilinear Higgs couplings, 
Xjihh an d Xhhh, and is the most useful one for extracting the coupling Xhhh- in Section 8 
we present, based on our calculations, the regions of the MSSM parameter space in which 
the trilinear couplings Xuhh and Xhhh could be measured; finally, in Section 9 we present 
a summary of our results and conclusions. 

2 The Higgs Sector of the MSSM 

In this section we review the Higgs sector of the Minimal Supersymmetric Standard Model 
in order to set the notation and to describe the approximations we use in our calculations. 
As mentioned in the introduction, we shall include the dependence on the parameters A 
and \x through mixing in the squark sector. Where there is an overlap, our notation and 
approach closely follow those of Ref . |lj . 

At the tree level, the Higgs sector of the MSSM is described by two parameters, 
which can be conveniently chosen as and tan/3 £Q. There are, however, substantial 
radiative corrections to the CP-even neutral Higgs masses and couplings :6]. In the one- 
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loop effective potential approximation, the radiatively corrected squared-mass matrix for 
the CP-even Higgs bosons can be written as [7j 



M' 



m\ sin 2 (3 + m| cos 2 (3 — (m| + m\) sin (3 cos (3 
— (m| + m\) sin (3 cos (3 m\ cos 2 (3 + m% sin 2 /3 



A n A 12 

A 12 A22 



(2.1: 



I 



where the second matrix represents the radiative corrections] 

The functions A^ depend, besides the top- and bottom-quark masses, on the Higgs 
bilinear parameter \x in the superpotential, the soft supersymmetry-breaking trilinear cou- 
plings (A t , A b ) and soft scalar masses (rriQ, my, m D ), as well as on tan/3. We shall ignore 
the 6-quark mass effects in in our calculations, which is a reasonable approximation 
for moderate values of tan (3 ^ 20-30. Furthermore, we shall assume, as is often done, 



A 



A t = A b , 



m = rrtQ = mjj = mo- 



(2.2) 



With these approximations we can write (m t is the top quark mass) [Zj: 



A 



11 



A 



22 



A 



12 



+ 



ml 



sin (3 
4 

sin 2 (3 
4 

sin 2 (3 



9 (ml, ml), 



fi(A + n cot (3) 

2 2 

m~ — ra- 
ti *2 



ml 



ra~ — ra~ raj 

tl 12 12 



+ fi cot /3) 

2 2 

m~ — ra- 
ti *2 



9 (ml, m 2 i2 ), 



m\ fi(A + fi cot /?) 
sin 2 /3 ra~ — ra~ 

' tl t2 



m f A(A + /icot/3) 2 2 

log — + 2 2 — 9(m- t ,m- 

ra~ mi — mi 1 2 

t2 ti t2 



(2.3) 



(2.4) 
(2.5) 



3 We note that two-loop corrections to the Higgs masses in the MSSM are sizable, especially for large 
mixing in the stop sector. For the dominant two-loop radiative corrections to the Higgs sector of the 
MSSM, see, e.g. Jjjj]. In this paper we restrict ourselves to one- loop corrections only. 
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where m~ and m~ 2 are squared stop masses given by 

™2 



mj = m t + m ± m t (A + fi cot f3) (2.6) 



*1,2 

(we have ignored the small D-term contributions to the stop masses) and 

m~ + m~ m~ 

tl C2 12 

The one-loop radiatively corrected masses (m/,, m#; < m#) of the CP-even Higgs 
bosons (h, H) can be obtained by diagonalizing the 2x2 mass matrix in Eq. (|2.1|) . The 
radiative corrections are, in general, positive, and they shift the mass of the lightest Higgs 
boson upwards from its tree-level value. We show in Fig. ^ the resulting mass of the 
lightest Higgs boson, m^, as a function of \i and tan /3, for two values of A and two values 
of 171a, and for fa — 1 TeV. With a wider range of parameter values, or when the squark 
mass scale is taken to be smaller, the dependence on \i and tan j3 can be more dramatic 
Q. 

The Higgs mass falls rapidly at small values of tan/3. Since the LEP experiments 
are obtaining lower bounds on the mass of the lightest Higgs boson, they are beginning 
to rule out significant parts of the small-tan (3 parameter space, depending on the model 



assumptions. For tan/5 > 1, ALEPH finds m h > 62.5 GeV at 95% C.L. [TUJ Q In 
our calculations, we shall therefore take tan (3 = 2 to be a representative value. [For a 
recent discussion on how the lower allowed value of tan j3 depends on some of the model 
parameters, see Ref. [T2"] .] 



Recently, a new study has been presented, with a lower limit of mh > 72.2 GeV, irrespective of tan (3, 
and a limit of - 88 GeV for 1< tan/3 < 2 [TT] . 
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3 Trilinear Higgs couplings 

The trilinear Higgs couplings that are of interest can be written jT3] as a sum of the 
tree-level coupling and one-loop radiative corrections: 



^Hhh 
^hhh 
^hAA 



^Hhh + A\ Hhh , 



hAA- 



(3.1) 
(3.2) 
(3.3) 



In units of gvnzj (2 cos 6*w) = (v / 2C j p) 1 / 2 m|, the tree-level couplings are given by 



x Hhh 



^1, 



hhh 



A hAA 



2 sin 2a sin(/3 + a) — cos 2a cos(/3 + a) 
3cos2asin(/5 + a), 
cos 2/3 sin(/3 + a), 



(3.4) 
(3.5) 
(3.6) 



with a the mixing angle in the CP-even Higgs sector, which can be calculated in terms of 
the parameters appearing in the CP-even Higgs mass matrix (j2.1j) . The one- loop radiative 
corrections in (j3.1j) - ()3.3j) are (in the above units): 



AA 



Hhh 



AA 



hhh 



3g 2 cos 2 9w m \ sin a cos 2 a 
16n 2 m^y sin 3 j3 



i i 

VTL r VTL r 

3 log -J^ + (ml 



III 



ml)C t {E t + 2F t )\og- 



VTLr 



ill 



H 



+2 3" [1 + « - OCtE t ] [I + (m| - ml 2 )C t F t }- 



+ 



2 
*2 



[1 



rrir — m 



3g 2 cos 2 9w TTT't cos3 a 
167T 2 m^, sin 3 /3 

2 2 2 

3 log + 3(m| - ml)C t F t log ' 1 



ml 



2 
*2 



(3.7) 



+ 2 ( ^ [1 + K - m|)CP 4 ] 3 + ™f- [l - « - m|)CP t ] 3 - 2 



3.8) 
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w 



sin 3 /? 



x log 



+ ( m |- m |)(^ + C t F t )log 



2 



+ (ml-mlfC t D 2 t F t g(ml,ml) 



(3.9) 



where 



D, 



E 



(A + n cot 13)/ '{ml -ml), 
(A - /xtan/3)/(m| - m|), 
(A + /icota)/(m| -m|), 
(A — n tan a) / (m| — m~ ) , 



(3.10) 



and we have ignored the contributions from 6-quarks and 6-squarks, which are in general 
small with respect to those arising from t-quarks and t-squarks. We have also adopted 
the simplification described in Eq. ()2.2|) in writing the above results. We shall make these 
approximations throughout this paper. 

We show in Figs. El El and |U the couplings Xnhh, Xhhh and XhAA as functions of /i 
and tan/?, for two values of A and two values of m^, all for m = 1 TeV. The explicit 
dependence on A and \i is not dramatic, but it should be kept in mind that unless is 
rather small, m^ may change considerably with A. 

The trilinear couplings change significantly with m^, and thus also with m^. This is 
shown more explicitly in Fig.0 where we compare Xnhh, Xhhh and Xhaa for three different 
values of tan/?, and the SM quartic coupling A SM . The SM quartic coupling includes one- 
loop radiative corrections |HJ, and its normalization is such that at the tree- level, it 
coincides with the trilinear coupling. 

At low values of m^, the MSSM trilinear couplings are rather small. For some value 
of m/j the couplings Xnhh and Xhhh start to increase in magnitude, whereas XhAA remains 
small. The values of m^ at which they start becoming significant depend crucially on 
tan (3. For tan (3 = 2 (Fig. |^) this transition takes place around m^ ~ 90-100 GeV, 
whereas for tan/3 = 5 and 15, the critical values of mh increase to 100-110 and 120 GeV, 
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respectively (see Figs. EJd and c). In this region, the actual values of \uhh and Xhhh (f° r a 
given value of rrih) change significantly if A becomes large and positive. A non-vanishing 
squark- mixing parameter A is thus seen to be quite important. Also, we note that for 
special values of the parameters, the couplings may vanish [T3]. See also Fig. 1 of Ref. [3]. 

To sum up the behaviour of the trilinear couplings, we note that Xuhh and Xhhh 
are small (< 1) for rrih ^ 100-120 GeV, depending on the value of tan (3. However, 
as rrih approaches its maximum value, which is reached rapidly as becomes large, 
rriA ^ 200 GeV, these trilinear couplings become large (~ 2 — 4). Thus, as functions of 
rriA, the trilinear couplings Xnhh an d Xhhh are large for most of the parameter space. We 
also note that, for large values of tan/5, Xnhh tends to be relatively small, whereas Xhhh 
becomes large, if also rriA (or, equivalently, m h ) is large. 

We note that for a given Higgs boson mass rrih, the tree level SM trilinear Higgs 
coupling is given by 

= 3(m h /m z ) 2 . (3.11) 

On the other hand, for large values of tjia (the decoupling limit) the corresponding MSSM 
trilinear coupling, Eq. ()3.5|) . becomes 

Xlhh = 3 cos(2a) sin(/3 + a) — > 3(m h /m z ) 2 , (3-12) 

i.e., it approaches the SM trilinear coupling. 

4 Production mechanisms 

The different mechanisms for the multiple production of the MSSM Higgs bosons in e + e~ 
collisions have been discussed by DHZ. The dominant mechanism for the production of 
multiple CP-even light Higgs bosons (h) is through the production of the heavy CP- 
even Higgs boson H, which then decays by H — > hh. The heavy Higgs boson H can 
be produced by P-strahlung, in association with A, and by the resonant WW fusion 
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mechanism. These mechanisms for multiple production of h 

e+e" -> ZH, AH 1 

\ , H ^hh, (4.1) 

e + e~ — > v e v e H j 

are shown in Fig. El We note that all the diagrams of Fig. El involve the trilinear coupling 

^Hhh- 

A background to (|4.1j) comes from the production of the pseudoscalar A in association 
with h and its subsequent decay to hZ 

e + e- ^hA, A — > ftZ, (4.2) 

leading to Zhh final states [see Fig. 0i] . 

A second mechanism for hh production is double Higgs-strahlung in the continuum 
with a Z boson in the final state [see Fig. E^-d], 

e + e~ -> Z* -> Z/i/i. (4.3) 

We note that the Feynman diagram of Fig.03 involves, apart from the coupling Xnhh, the 
trilinear Higgs coupling as well, whereas the other diagrams do not involve any of 
the trilinear Higgs couplings. 

A third way of generating multiple Higgs bosons in e + e~ collisions is through associated 
production of (hh) with the pseudoscalar A in the continuum [see Fig. [BJ: 

e + e" -> -> /i/iA. (4.4) 

This process will be briefly discussed in Section 6. It involves, besides Xnhh and Xhhh, the 
trilinear coupling XhAA as well. It is, however, difficult [4 a to measure this coupling XhAA 
through the process ()4.4j) . 

Finally, there is a mechanism of multiple production of the lightest Higgs boson 
through non-resonant WW (ZZ) fusion in the continuum [see Fig. [Hj: 

e+e" -> z/^eiy-W* -> P e i/ e /i/i, (4.5) 

which will be discussed in Section 7. 
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It is important to note that all the diagrams of Fig. |H1 involve the trilinear coupling 
\jihh only. On the other hand, Fig. [7^, Fig. |HJd and Fig. all involve both the trilinear 
Higgs couplings X Hhh and X hhh . 



by the cascade decay H — > hh. In terms of the Z-electron couplings v e = —1 + 4 sin 2 6 



5 Higgs-strahlung and Associated Production of H 

As stated in Section 4, the dominant source for the production of multiple Higgs bosons 
(h) in e + e~ collisions is through the production of the heavier CP-even Higgs boson H 
either via Higgs-strahlung or in association with A |T|, followed, if kinematically allowed, 

n- 

a e = —1, the cross sections for these processes can be written as [T7)| IT7] 

, i „„ s G%mt , 9 9 . „.„ , AV 2 \\y + 12m 2 7 /s] 

where Aj refers to A(m 2 ,m^; s), the two-body phase-space function, and is given as 



AK,m;m)= 1-^-4 - (5.3) 

In Fig. EI we plot the cross sections (jQJ) and ()5.2)1 for the e + e~ centre-of-mass energies 
^/i = 500 GeV and 1.5 TeV, as functions of the Higgs mass m# and for tan/? = 2.0. For 
large values of the mass m^ of the pseudoscalar Higgs boson, all the Higgs bosons, except 
the lightest one (h), become heavy and decouple [TH] from the rest of the spectrum. In 
this case 

i, n x mi sin 2 4/? 
cos 2 03 - a) ~ ^ ; (5.4) 

and the associated AH production (J5.2)) becomes the dominant production mechanism 
for H. 

At values of tan/? that are not too large, the trilinear Hhh coupling Xuhh can be 
measured by the decay process H — ► ftA, which has a width 

G F>? Hhh mi ( _ Ami x 1 ' ? 



T(H^hh) = — H ™ M l ^ . (5.5) 

\Qtx^2 m H \ m 2 H ' 
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However, this is possible only if the decay is kinematically allowed, and the branching ratio 
is sizeable. In Fig. we show the branching ratios (at tan (3 = 2) for the main decay 
modes of the heavy CP-even Higgs boson as a function of the H mass. Apart from the 
hh decay mode, the other important decay modes are H —> WW*, ZZ*. (We have here 
disregarded decays to supersymmetric particles: charginos, stops, etc. If such particles 
are kinematically accessible, the H — > hh and A — > Zh rates could be much smaller [19].) 
We note that the couplings of H to gauge bosons can be measured through the production 
cross sections for e + e~ — > v e v e H\ therefore the branching ratio BR(H — > hh) can be used 
to measure the triple Higgs coupling Xnhh- 

The Higgs-strahlung process [Fig. |HK, Eq. 1)5.1)1 ] gives rise to resonant two-Higgs [hh] 
final states. This is to be contrasted with the associated production process [Fig. Eb; 
Eq. 1)5.2)1 ]. which typically yields three Higgs h[hh] final states, since the channel A — ► /iZ 
is the dominant decay mode of A in the mass range of interest. The decay width for 
A — > /iZ can be written as (2Hj 

T(A - = cos 2 (/5 - «)^A 1 / 2 ( m |, m * ; m 2 ) A(m 2 m\- m|), (5.6) 

87rV2 "^a 

where the A are phase-space factors given by Eq. ([5.3)1 . In Fig. ^2 we show the branching 
ratios for the pseudoscalar A for tan/5 = 2.0. 

A background to the multiple production of lighter Higgs bosons h comes from hh 
states generated in the sequential reaction e + e~ —>■ Ah — > [Zh]h [see Fig. [)Jl]. This is 
a genuine background in the sense that no Higgs self-couplings are involved. But these 
background events are expected to be topologically very different from the signal events, 
since the two h bosons do not form a resonance, whereas the [Zh] does. The cross section 
for the process e + e~ Ah can be written as [1.7] 

a(e+e- - Ah) = ^(, e 2 + aj) cos 2 (/? - a /^™^ , (5.7) 
9o7rs (1 — m z /s) z 

and is shown in Fig. Altogether with the signal cross sections ()5.1)) and (|5.2|) . As a 
consequence of the decoupling theorem ^H], the cross section becomes small for large 
values of ma- 
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For increasing values of tan /3, the Hhh coupling gradually gets weaker (see Figs. |2]and 
EJ), and hence the prospects for measuring Xuhh diminish. This is indicated by Fig. [T2l 
where we show the H and A branching ratios for tan/3 = 5. 

There is in fact a sizeable region in the m^-tan/3 plane where the decay H —>■ hh 
is kinematically forbidden. This is indicated in Fig. ED in this figure we also display 
the regions where the H — > hh branching ratio is in the range 0.1-0.9. Clearly, in the 
forbidden region, the Xnhh cannot be determined from resonant production. 

6 Double Higgs-strahlung and Triple h Production 

For small and moderate values of tan/3, the study of decays of the heavy CP-even Higgs 
boson H provides a means of determining the triple-Higgs coupling Xnhh- In order to 
extract the coupling Xhhh, other processes involving two-Higgs (h) final states must be 
considered. The Zhh final states, which can be produced in the double Higgs-strahlung 
e + e~ — ► Zhh of Fig. could provide one possible opportunity, since it involves the 
coupling Xhhh through the mechanism of Fig. 03. In this section we shall study these 
non-resonant processes in detail. 

6.1 The Double Higgs-strahlung e + e _ — > Zhh 

The doubly differential cross section for the process e + e~ — ► Zhh shown in Fig. [7|can be 
written as j~D 

da(e + e- -> Zhh) G%m% , 2 2 A 

dx x dx 2 384x/2vr3 s [Ve + aJ (1 - /iz) 2 ' 1 } 

where the couplings v e and a e have been defined at the beginning of Section 5. Because of 

some misprints in the formulas given in |3] for the coefficient A, we have recalculated it. 

Following jl], we introduce = 2Ei i2 / 'y/s for the scaled energies of the Higgs particles, 

x 3 = 2 — x\ — x 2 for the scaled energy of the Z boson, and = 1 — x^. Also, we denote 

by ^ = mf/s the scaled squared masses of various particles: 

fih = m 2 h /s, fiH = m 2 H /s, Hw = m 2 w /s. (6.2) 
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We can express our result in a compact form as follows: 



A = fi Z {±M 2 f a + |%i)| 2 f b + 2Re[ab*( yi )} g ab + Re[b( Vl )b* (y 2 )} g bb ] + { Xl «-> x 2 }. 

(6.3) 



Here, 



sin(/3 - a)\ hhh cosjfi - a)X H hh 

yz + hz - P>h ys + nz- he 



+ 



sin 2 (/? — a) sin 2 (/? — a) 



yx + hk - ^z V2 + Hh - P-z 



+ 



(6.4) 



represents a contribution from diagram [T^i, where the lepton tensor couples directly to 
the final-state Z polarization tensor, as well as the contributions of diagrams [7)d andEk- 
Similarly, 



Ky) 



1 ( sin 2 (/5 — a) cos 2 (/3 — a) 



(6.5) 



2/xz \y + Hh ~ Hz y + Hh- Ha, 

represents the part of diagram UK where the lepton tensor couples to the final-state Z 
polarization tensor indirectly via the Higgs momenta q± and q 2 , as well as diagram [7|i. 
The tildes on keep track of the widths, e.g. Hz = (w§ + im z T z )/s. 

The Higgs self- couplings Xuhh and Xhhh occur only in the function a, Eq. ()6.4|) . The 
coefficients / and g>, which do not involve any Higgs couplings, can be expressed rather 
compactly as 



f a = x 2 + 8fl Z , 

fb = (xl - A/j, h ) [(y x - Hz) 2 - ^HzHh] , 

g a b = Hz[2(hz - ^Hh) + x\ + x 2 (x 2 + x 3 )} - yi(2y 2 - x x x 3 ), 
9bb = Hzi^Hh + 6 - xix 2 ) + 2\x z (h\ + 2/3 - 4/i/J 

+(j/3 - x±x 2 - x 3 fx z - ^HhHz)(2y 3 - x x x 2 - 4^ h + 4/i Z ) 



(6.6) 



These coefficients (we use a mixed notation, which involves both x and y) correspond to 
those of |lj as follows: (f a , f b ,g a b,gbb) = (fo, fx, fs, $2)- With this identification, we agree 
with the result given in the Erratum to [I]. 
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In the limit of large uia, sin(/3 — a) — > 1, the cross section reduces to the Standard 
Model cross section with 



^hhh 



2 2/3 + ^ - P>h 



1 1 

+ 



J/i + - fiz V2 + Hh - 



+ Wz (6 ' 7) 



%) = ^~ 7 1 (6-8) 

where at the tree level, Xhhh —> X^ h , as discussed in Sec. 3. 

We show in Fig. Hlk the Z/i/i cross section, as given by Eqs. (|5.1j) . (|5.5jl and (|6.1jl . in 
the limit of no squark mixing, and with m = 1 TeV. The structure around = 70 GeV 
is due to the vanishing and near- vanishing of the trilinear coupling. In Fig. ITlb-d we have 
introduced squark mixing: A = 1 TeV, \x = 0, ±1 TeV. (For the decoupling-limit cross 
section, which is also shown, we use the MSSM coupling, instead of the SM coupling, for 
the reason given in Sec. 3.) 

In the case of no mixing, there is a broad minimum from rrih — 78 to 90 GeV, followed 
by an enhancement around ~ 90-100 GeV. This structure is due to the vanishing of 
the branching ratio for H — > hh, which is kinematically forbidden in the region ~ 78- 
90 GeV, see Fig. El (this coincides with the opening up of the channel H — > WW), 
followed by an increase of the trilinear couplings. This particular structure depends 
considerably on the exact mass values m# and m^. Thus, it depends on details of the 
radiative corrections and the mixing parameters A and \x. 

The A — ► hZ channel contributes of the order of 20% in the region of the maximum 
at m h = 90-100 GeV. 



6.2 Triple-/^ production 

The resonant and non-resonant production of Ahh [Fig. |H] can lead to three-/;, final states 
in the region of tua, where A has a significant branching ratio for decaying to Zh, i.e. for 
rriA below the ti threshold, and relatively low values of tan j3 [cf . Figs. and E] • 

In principle, this channel allows for a study of the coupling Xhaa [cf • Fig. [S^,] . However, 
the prospects for measuring this coupling, which is rather small [see Fig. Hj, was studied 
in Ref. |lj and found not to be very encouraging. 
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7 Fusion Mechanism for Multiple-/^ Production 

As mentioned in Section 4, a double Higgs (hh) final state in e + e~ collisions can also 
result from the WW fusion mechanism, which can either be a resonant process as in 
(I4.1J1 . or a non-resonant one like (|4.5|) . Since the neutral-current couplings are smaller 
than the charged- current ones, the cross section for the ZZ fusion mechanism in and 
(14. 5j) is an order of magnitude smaller than the WW fusion mechanism. We shall thus, 
in the following, ignore the ZZ fusion mechanism, and concentrate instead on the WW 
mechanism. 

7.1 Resonant WW fusion 



The WW fusion mechanism provides another large cross section for the multiple produc- 
tion of h bosons. The cross section for e + e~ — > Hv e v e can be written as J2T] 



cr(e + e~ 
where 



Hv e v e ) 



64V2 



7P 



dx 



dy 



x [l + (y-x)/nwY 



F{x,y) 



cos 2 (/?-a), (7.1) 



F(x,y) = 16[F{x,y) + G{x,y)} 
F(x,y) = 



G{x,y) 



2x l + 2x 2 + x 1 

y 3 y 2 2y 2 



x 2 + x 1 

n + 



v 2y 2 

with fii defined by Eq. (|fj.2jl and 



l + z 



l + z 
log(l + z 



logfl + z) 



yjx - n H ) 



(7.2) 

f4 ^ (1 ~ y) ; (7.3) 



2/ 3 (l + z) 



(7.4) 



(7.5) 



For -y/i, mn ^> m^, and in the effective longitudinal W approximation, the cross section 
(17. lj) for e + e~ — > Hv e v e can be written in the following simple form [2*2*] 



trie e 



G 3 F m^ 



AV2 



7T" 



1 + 1 Iok — - 2 [ 1 



m 



H 



m 



li 



cos 2 (/5-a). (7.6) 



However, in this approximation the cross section may be overestimated by a factor of 2 
for small values of masses and/or small centre-of-mass energies. For example, at \fs = 



16 



500 GeV the equivalent W approximation gives a result that is twice as large as the exact 
cross section. Therefore, we use the exact cross section (j7.1j) in our calculations. 

The cross section (17. lj) is plotted in Fig. EH for centre-of-mass energies, yfs = 500 GeV 
and 1.5 TeV, and for tan/5 = 2.0, as a function of tuh- The resonant fusion mechanism, 
which leads to [hh] + [missing energy] final states is competitive with the process e + e~ — > 
HZ — > [hh] + [missing energy], particularly at high energies. Since the dominant decay 
of h will be into bb pairs, the if-strahlung and the fusion mechanism will give rise to 
final states that will predominantly include four 6-quarks. On the other hand, the process 
e + e _ — > AH will give rise to six 6-quarks in the final state, since the AH final state 
typically yields three-Higgs h[hh] final states. 

7.2 Non-resonant fusion WW — > hh 

Besides the resonant WW fusion mechanism for the multiple production of h bosons, 
there is also a non-resonant WW fusion mechanism: 

e + e~ — > v e i> e hh, (7.7) 

through which the same final state of two h bosons can be produced. The cross section 
for this process, which arises through WW exchange as indicated in Fig. El can be written 
in the 'effective WW approximation' a|] 



dx — a ww (x), (7.8) 

where r = 4m\/s. In the above, the cross section is written as a WW cross section, at 
invariant energy squared s = xs, folded with the WW 'luminosity' [22J: 

dL{x) G\m A w fv 2 + a 2 \ 2 1 f 1 



da; 2 \ An 2 J x \ x 

where v 2 + a 2 = 2. 



-<m +x)log--2(l-x) , (7.9) 



3 There could be sizable deviations of the effective WW approximation from the exact result. 
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The WW cross section receives contributions from several amplitudes, according to 
the diagrams (a)-(d) express the 

result in a form analogous to that of Ref. 



2 



fi z sm(/3-a) [i z cos(f3 - a) 

: ; Ahhh H : ; *Hhh + ± 

1- Hh 1 - 

fi z sm(p-a) fi z cos(/3 - a) ' 

: ; Ahhh H ; ; AHhh + 1 

1 - n h 1 - fi H 

x [sin 2 (/? — a) g\ + cos 2 (/3 — a) #2] 



+ {sin 4 (/? - a) 33 + cos 4 (/5 - a) 34 + sin 2 [2(/3 - «)] s 5 } L (7.10) 
where we have introduced 'reduced squared masses' 

fi z = m 2 z /s, fi w =m 2 w /s, fl h = m 2 h /s, fi H = m 2 H /s, (7.11) 

1 /2 

and the Higgs velocity is (3h — (1 — /i^) 

Our approach differs from that of DHZ in that we do not project out the longitudinal 
degrees of freedom of the intermediate W bosons. Instead, we follow the approach of 
Ref. 23J, where transverse momenta are ignored everywhere except in the W propagators, 
the integrations over which are approximated as (here p\ and p[ denote electron and 
neutrino momenta, respectively, in the process (J7.7j) ): 

d 2 ^ n ^ ~ X ^ (7 12) 

Pl± [(pi -Pi) 2 ~ m w} 2 ~ m w 

where Xi [and x 2 ] represents the energy of the W. The energy squared of the subprocess 

is given as s = (p\Xi + p 2 x 2 ) 2 = Xix 2 s = xs. 

6 For each of the diagrams (a) and (d) there are two contributions, corresponding to the interchange 
of the two Higgs particles. 

7 There are some misprints in Ref. 0], so we present here results of an independent calculation. 

8 In Ref. 0], the factor in front of the term involving g\ and 172 reads fiw/Wwfih)] it should be 
(1 + fiw)/(PwPh)- With (3w — 1, the prefactors in their Eq. (16) reduce to ours. 
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The contributions of diagrams (b) + (c), (a) and (d) are given by the terms go, g 3 and 
g 4 , respectively, with go(x) = 1, and 

g 3 (x) = 8/3 h [2p, w + (fi h - jj, w f}[2fL W + 1 - 3(p, h - fiwf] — 
+16[2/V + (P-h ~ frwYYvw + 16$j(l + a w ) 2 , 

g 4 (x) = 8f3 h (p, h - /ic) 2 [l - 3(// fc - Ac) 2 ] — 

+16(/tfc - Ac) 4 yc + 16^(1 + a c ) 2 , (7.13) 

where 

/W = l0g l- 2/ 4 + 2^ + ^ (7 ' 14) 



(1 - 2^ + 2Aw) 2 



= ^ o., ,Z. ^_ R 2 i ( 7 - 15 ) 



and similarly l c , yc and ac, with \x w replaced by jlc, the latter being defined in terms of 
the charged Higgs mass m H +. 

The interference between diagrams (b) + (c) and (a) is given by the term gx, whereas 
the interferences between diagrams (b) + (c) and (d), and between (a) and (d) are given 
by g% and g$, respectively. For these interference terms, we find 

gi(x) = 8[2fi w + {fth - frw) 2 ]lw - 4/3^(1 + 2fx h - 2p, w ), 
g 2 {x) = 8(jl h -jlc) 2 lc-4:P h (l + 2jl h -2jl c ), 

9b(x) = ^(Z w lw + Z c lc) + 8(3 2 h (l + a w )(l + a c ), (7.17) 

with 

Z w = (1 + 2aw)2 [8fi w + (1 + 2a w f] + (1 ~ 2aw)2 [SfM W + (1 + 2a w )% (7.18) 

CLc ~ Ci\V a C + a w 

(l + 2a c ) 2 r , s9l (l + 2a c ) 2 r , s9l , 

Z c = ~- —Ww + (1 + 2ac) 2 ] + ^ ~^[8jl w + (l-2ac) 2 }. (7.19) 

CLc — a w CLc T CLw 

Our functions g\-g<z> correspond to those of DHZ, cf. our Eq. ()7.1()|) and their Eq. (16) 
in jlj. At small m h , the cross section is sensitive to small x, where the 'effective WW 
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approximation' is not well defined, and our results differ from those of DHZ. However, 
apart from the contributions from small x, our results agree with those of DHZ to a 
precision of 1-5%. 

We show in Fig. [T3] the WW fusion cross section, at y/s = 1.5 TeV, as given by 
Eqs. f|7.1|) and (|7.8|) . in the limit of no squark mixing, as well as with mixing (as indicated), 
and with rh = 1 TeV. The structure is very reminiscent of that of Fig. El and for the 
same reasons. However, the scale is different. 

8 Sensitivity to X Hhh and X hhh 

Following 4J, we have indicated in the m^-tan/3 plane the regions where Xnhh and Xhhh 
might be measurable, according to criteria analogous to those given there. In Fig. El we 
consider y/s = 500 GeV, and identify regions according to the following criteria: 

(i) Regions where XHhh might become measurable are identified as those where cr(H) x 
BK(H — > hh) > 0.1 fb (solid), with the simultaneous requirement of 0.1 < BK(H — > 
hh) < 0.9 [see Figs. ITTVll3j . In view of the recent, more optimistic, view on the 
luminosity that might become available, we also give the corresponding contours for 
0.05 fb (dashed) and 0.01 fb (dotted). For a(H) we take the sum of (jSHJ), (Q and 
(EH). 

(ii) Regions where Xhhh might become measurable are those where the continuum WW — 
hh cross section [Eq. ()7.8|) ] is larger than 0.1 fb (solid). Also included are contours 
at 0.05 (dashed) and 0.01 fb (dotted). 

Such regions are given for four cases of the mixing parameters A and fi, as indicated. We 
have excluded from the plots the region where nih < 62.5 GeV, according to the LEP 
lower bound |10j . This corresponds to low values of m^. 

We note that with an integrated luminosity of 500 fb -1 , the contours at 0.1 fb cor- 
respond to 50 events per year. This will of course be reduced by efficiencies, but should 
indicate the order of magnitude that can be reached. 
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At \/s = 500 GeV, with a luminosity of 500 fb -1 per year, the trilinear coupling Xuhh 
is accessible in a considerable part of the m^-tan/3 parameter space: at of the order 
of 200-300 GeV and tan/3 up to the order of 5. With increasing luminosity, the region 
extends somewhat to higher values of m^. 

At values of tua below 100 GeV, there is also a narrow band where Xuhh is accessible. 

The 'steep' edge around uia — 200 GeV (where increased luminosity does not help) is 
determined by the vanishing of BR(H — > hh), see Fig. fTTTl 

The coupling Xhhh is accessible in a much larger part of this parameter space, but with 
a moderate luminosity, 'large' values of tan (3 are accessible only if A is small. 

In Fig. H3 we consider y/s = 1.5 TeV, and present the analogous contours. Here, for 
the case of Xnhh we demand cr(H) x BR(iJ — > hh) > 0.5 fb (solid) and 0.1 fb (dashed), and 
for the case of Xhhh we require the corresponding cross section [Eq. ([7.8)1 ] to be larger than 
0.5 fb (solid) and 0.1 fb (dashed). If a luminosity corresponding to these cross sections 
becomes available at ^/s = 1.5 TeV, a somewhat larger region than at ^/s = 500 GeV is 
accessible in the m^-tan (3 plane. 

It should be stressed that the requirements discussed here are necessary, but not suf- 
ficient conditions for the trilinear couplings to be measurable. We also note that there 
might be sizable corrections to the WW approximation, and that it would be desirable 
to incorporate the dominant two-loop corrections to the trilinear couplings in the calcu- 
lations. 

9 Conclusions 

We have carried out a detailed investigation of the possibility of measuring the MSSM 
trilinear couplings Xrhh and Xhhh at an e + e~ collider. Where there is an overlap, we have 
confirmed the results of Ref. j3j. Our emphasis has been on taking into account all the 
parameters of the MSSM Higgs sector. We have studied the importance of mixing in the 
squark sector, as induced by the trilinear coupling A and the bilinear coupling \i. 

At moderate energies (y/s = 500 GeV) the range in the m^-tan (3 plane that is accessi- 
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ble for studying Xuhh changes quantitatively for non-zero values of the parameters A and 
/i. As far as the coupling \ hhh is concerned, however, there is a qualitative change from 
the case of no mixing in the squark sector. If A is large, then high luminosity is required to 
reach 'high' values of tan/3. At higher energies (y/s = 1.5 TeV), the mixing parameters A 
and \i change the accessible region of the parameter space only in a quantitative manner. 
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Figure 1: Mass of the lightest Higgs boson mt as a function of fi and tan/5. Two values 
of vriA and two values of A are considered: a) itla = 100 GeV, A — 0, b) itla = 400 GeV, 
A = GeV, c) m A = 100 GeV, A = 1 TeV, d) m A = 400 GeV, A = 1 TeV. We have 
taken fh = 1 TeV. 
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Figure 2: Trilinear Higgs coupling Xuhh as a function of \i and tan/3. The values of the 
parameters are the same as in Fig. ^ 
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Figure 3: Trilinear Higgs coupling X hhh as a function of fi and tan/?. The values of the 
parameters are the same as in Fig. ^ 
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Figure 4: Trilinear Higgs coupling XhAA as a function of \x and tan (5. The values of the 
parameters are the same as in Fig. ^ 
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Figure 5: Trilinear Higgs couplings \hhh-i ^hhh and Xhaa as functions of m h for three 
values of tan/?: (a) tan/? = 2.0, (b) tan/? = 5.0, (c) tan/? = 15. Each coupling is shown 
for three cases of the mixing parameters: no mixing (A = 0, fi = 0, solid), mixing with 
A = 1 TeV and fi = -1 TeV (dotted), as well as A = 1 TeV and n = 1 TeV (dashed). 
For comparison the SM quartic coupling A SM is also shown. 
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Figure 6: Feynman diagrams for the resonant production of hh final states in e + e~ colli- 
sions. Diagrams (a) and (b) represent the production of H in association with Z and A, 
respectively, whereas diagram (c) is the WW fusion mechanism for the production of H. 
The Higgs boson H decays via H — > hh to produce the two-Higgs final state. 




Figure 7: Feynman diagrams for the non-resonant production of hh final states in associ- 
ation with Z. The diagram (d), with A produced on the mass shell, which subsequently 
decays via A — > hZ, is a background to the resonance process, Fig. 




Figure 8: Feynman diagrams for the associated production of hh with the pseudoscalar 
A in the continuum. 




Figure 9: Feynman diagrams for the non-resonant WW fusion mechanism for the pro- 
duction of hh states in e + e~ collisions. 
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Figure 10: Cross sections for the production of the heavy Higgs boson H in e + e" collisions. 
Also shown is the cross section for the background process in which Ah is produced in the 
final state. We have taken y/s = 500 GeV and 1.5 TeV. Solid curves are for no mixing, 
A — 0, /j, — 0. Dashed and dotted curves refer to mixing: A = 1.0 TeV, /j, — 1.0 TeV 
(dashed) and A = 1.0 TeV, yt, = -1.0 TeV (dotted). 
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Figure 11: Branching ratios for the decay modes of the CP-even heavy Higgs boson H, 
and the CP-odd Higgs boson A for tan (3 = 2.0. Solid curves are for no mixing, A = 0, 
fi = 0. Dashed and dotted curves refer to mixing: A = 1.0 TeV, /i = 1.0 TeV (dashed) 
and A = 1.0 TeV, fi = -1.0 TeV (dotted). 
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Figure 12: As in Fig. ^2 but with tan/3 = 5.0. Solid curves are for no mixing, A = 0, 
\i = 0. Dashed and dotted curves refer to mixing: A = 1.0 TeV, /i = 1.0 TeV (dashed) 
and A = 1.0 TeV, /i = -1.0 TeV (dotted). 
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Figure 13: The region in the m^-tan/3 plane where the decay H — *■ hh is kinematically 
forbidden is indicated by a solid line contour. Also given are contours at which the 
branching ratio equals 0.1 (dotted), 0.5 (dashed) and 0.9 (dash-dotted, at the far left). 
Four cases of mixing parameters A and \x are considered, as indicated. 
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Figure 14: Cross section a(e + e~ — > Zhh) as a function of m h) for four cases: (a) no 
mixing, (b)-(d) A = 1 TeV, /x = 0, —1 and 1 TeV, as indicated. The dotted curve is the 
resonant production. The dashed curve gives the decoupling limit. 
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Figure 15: Cross section a(e + e~ — > v e i) e hh) (via WW fusion) as a function of m h , for 
four cases: (a) no mixing, (b)-(d) A = 1 TeV, /i — 0, — 1 and 1 TeV, as indicated. The 
dotted curve is the resonant production. The dashed curve gives the decoupling limit. 
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Figure 16: Regions where trilinear couplings Xuhh an d Xhhh might be measurable at 
yfe = 500 GeV. Inside contours labelled X Hhh , a{H) x BR(H — > hh) > 0.1 fb (solid), 
while 0.1 < BR(if — > /i/i) < 0.9. Inside (to the right or below) contour labelled Xhhh, the 
continuum WW — > fo/i cross section exceeds 0.1 fb (solid). Analogous contours are given 
for 0.05 (dashed) and 0.01 fb (dotted). Four cases of mixing are considered, with A = 
or 1 TeV, and fi — or ±1 TeV, as indicated. 
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Figure 17: Regions where trilinear couplings Xuhh an d Xhhh might be measurable at 
s/s = 1.5 TeV. The contours correspond to 0.5 fb (solid) and 0.1 fb (dashed). All other 
parameters are the same as in Fig. ED 
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